A comparative study on the detection of Mycobacterium leprae DNA in urine samples of leprosy patients using Rlep-PCR with other conventional samples.

BACKGROUND: Mycobacterium leprae causes leprosy that is highly stigmatized and chronic infectious skin disease. Only some diagnostic tools are being used for the identification M. leprae in clinical samples, such as bacillary detection, and histopathological tests. These methods are invasive and often have low sensitivity. Currently, the PCR technique has been used as an effective tool fordetecting M. leprae DNA across different clinical samples. The current study aims to detect M. leprae DNA in urine samples of untreated and treated leprosy patients using the Rlep gene (129 bp) and compared the detection among Ridley-Jopling Classification. METHODS: Clinical samples (Blood, Urine, and Slit Skin Smears (SSS)) were collected from leprosy and Non-leprosy patients. DNA extraction was performed using standard laboratory protocol and Conventional PCR was carried out for all samples using Rlep gene target and the amplicons of urine samples were sequenced by Sanger sequencing to confirm the Rlep gene target. RESULTS: The M. leprae DNA was successfully detected in all clinical samples across all types of leprosy among all the study groups using RLEP-PCR. Rlep gene target was able to detect the presence of M. leprae DNA in 79.17% of urine, 58.33% of blood, and 50% of SSS samples of untreated Smear-Negative leprosy patients. The statistical significant difference (p = 0.004) was observed between BI Negative (Slit Skin Smear test) and RLEP PCR positivity in urine samples of untreated leprosy group. CONCLUSION: The PCR positivity using Rlep gene target (129 bp) was highest in all clinical samples among the study groups, across all types of leprosy. Untreated tuberculoid and PNL leprosy patients showed the highest PCR positivity in urine samples, indicating its potential as a non-invasive diagnostic tool for leprosy and even for contact screening.

Mycobacterium leprae and host immune transcriptomic signatures for reactional states in leprosy.

Background: Mycobacterium leprae transcriptomic and human host immune gene expression signatures that demonstrate a plausible association with type I (T1R) and type II reactions (T2R) aid in early diagnosis, prevention of nerve damage and consequent demyelinating neuropathy in leprosy. The aim of the study is to identify M. leprae and host-associated gene-expression signatures that are associated with reactional states in leprosy. Methods: The differentially expressed genes from the whole transcriptome of M. leprae were determined using genome-wide hybridization arrays with RNA extracted from skin biopsies of 20 T1R, 20 T2R and 20 non reactional controls (NR). Additionally, human immune gene-expressions were profiled using RT2-PCR profiler arrays and real-time qPCRs. Results: The RNA quality was optimal in 16 NR, 18 T1R and 19 T2R samples. Whole transcriptome expression array of these samples revealed significant upregulation of the genes that encode integral and intrinsic membrane proteins, hydrolases and oxidoreductases. In T1R lesional skin biopsy specimens, the top 10 significantly upregulated genes are ML2064, ML1271, ML1960, ML1220, ML2498, ML1996, ML2388, ML0429, ML2030 and ML0224 in comparison to NR. In T2R, genes ML2498, ML1526, ML0394, ML1960, ML2388, ML0429, ML0281, ML1847, ML1618 and ML1271 were significantly upregulated. We noted ML2664 was significantly upregulated in T1R and repressed in T2R. Conversely, we have not noted any genes upregulated in T2R and repressed in T1R. In both T1R and T2R, ML2388 was significantly upregulated. This gene encodes a probable membrane protein and epitope prediction using Bepipred-2.0 revealed a distinct B-cell epitope. Overexpression of ML2388 was noted consistently across the reaction samples. From the host immune gene expression profiles, genes for CXCL9, CXCL10, CXCL2, CD40LG, IL17A and CXCL11 were upregulated in T1R when compared to the NR. In T2R, CXCL10, CXCL11, CXCL9, CXCL2 and CD40LG were upregulated when compared to the NR group. Conclusion: A gene set signature involving bacterial genes ML2388, ML2664, and host immune genes CXCL10 and IL-17A can be transcriptomic markers for reactional states in leprosy.

Molecular epidemiology and transmission dynamics of leprosy among multicase families and case-contact pairs.

BACKGROUND: Human-to-human transmission of Mycobacterium leprae among household contacts of active leprosy cases is significant, and surveillance of household contacts is vital to interrupting the transmission chain for this disease. This study was conducted to identify similarities in M. leprae strains, based on genomic single nucleotide polymorphisms (SNPs), among cases and their household contacts and in multicase families in order to decipher possible associations, transmission links, various clinical conditions of index cases that enhance person-to-person transmission, and timelines for transmission patterns. METHODS: PCR for M. leprae DNA detection (amplification of the Rlep gene) and SNP subtyping of M. leprae strains was performed for 61 index cases and one of their household contacts. Additionally, we studied six families with multiple cases of leprosy, to understand timelines of infectivity and its relation to severity of the disease in the index cases. RESULTS: Index cases with lepromatous (LL) and borderline lepromatous (BL) leprosy, together with a positive bacteriological index (BI) for M. leprae, result in a higher percentage of their contacts subclinically infected with M. leprae, with odds ratios (OR) of 6.6 (95% confidence interval (CI) 1.6-27.6) for BL and LL, and 7.07 (CI 1.41-35.41) for BI-positive index cases. 75% of the case-contact pairs had a similar SNP subtype of M. leprae. The timeline of infection in multicase families revealed that contacts were infected during the BI-positive period of the index case. CONCLUSION: Using molecular methods, we determined that positivity for M. leprae DNA in contacts of index leprosy cases was attributed to clinical characteristics of leprosy in the index cases. LL and BL forms of leprosy, together with positive BI, contributed to dissemination of infection to household contacts. In conclusion, we found a relationship between SNP subtypes within index case-contact pairs. This method can help decipher the transmission patterns and identify individuals at risk of contracting leprosy.

Computer-aided synthesis of dapsone-phytochemical conjugates against dapsone-resistant Mycobacterium leprae.

Leprosy continues to be the belligerent public health hazard for the causation of high disability and eventual morbidity cases with stable prevalence rates, even with treatment by the on-going multidrug therapy (MDT). Today, dapsone (DDS) resistance has led to fear of leprosy in more unfortunate people of certain developing countries. Herein, DDS was chemically conjugated with five phytochemicals independently as dapsone-phytochemical conjugates (DPCs) based on azo-coupling reaction. Possible biological activities were verified with computational chemistry and quantum mechanics by molecular dynamics simulation program before chemical synthesis and spectral characterizations viz., proton-HNMR, FTIR, UV and LC-MS. The in vivo antileprosy activity was monitored using the 'mouse-foot-pad propagation method', with WHO recommended concentration 0.01% mg/kg each DPC for 12 weeks, and the host-toxicity testing of the active DPC4 was seen in cultured-human-lymphocytes in vitro. One-log bacilli cells in DDS-resistant infected mice footpads decreased by the DPC4, and no bacilli were found in the DDS-sensitive mice hind pads. Additionally, the in vitro host toxicity study also confirmed that the DCP4 up to 5,000 mg/L level was safety for oral administration, since a minor number of dead cells were found in red color under a fluorescent microscope. Several advanced bioinformatics tools could help locate the potential chemical entity, thereby reducing the time and resources required for in vitro and in vitro tests. DPC4 could be used in place of DDS in MDT, evidenced from in vivo antileprosy activity and in vitro host toxicity study.

Computational saturation mutagenesis to predict structural consequences of systematic mutations in the beta subunit of RNA polymerase in Mycobacterium leprae.

Rifampin resistance in leprosy may remain undetected due to the lack of rapid and effective diagnostic tools. A quick and reliable method is essential to determine the impacts of emerging detrimental mutations in the drug targets. The functional consequences of missense mutations in the ß-subunit of RNA polymerase (RNAP) in Mycobacterium leprae (M. leprae) contribute to phenotypic resistance to rifampin in leprosy. Here, we report in-silico saturation mutagenesis of all residues in the ß-subunit of RNAP to all other 19 amino acid types (generating 21,394 mutations for 1126 residues) and predict their impacts on overall thermodynamic stability, on interactions at subunit interfaces, and on ß-subunit-RNA and rifampin affinities (only for the rifampin binding site) using state-of-the-art structure, sequence and normal mode analysis-based methods. Mutations in the conserved residues that line the active-site cleft show largely destabilizing effects, resulting in increased relative solvent accessibility and a concomitant decrease in residue-depth (the extent to which a residue is buried in the protein structure space) of the mutant residues. The mutations at residue positions S437, G459, H451, P489, K884 and H1035 are identified as extremely detrimental as they induce highly destabilizing effects on the overall protein stability, and nucleic acid and rifampin affinities. Destabilizing effects were predicted for all the clinically/experimentally identified rifampin-resistant mutations in M. leprae indicating that this model can be used as a surveillance tool to monitor emerging detrimental mutations that destabilise RNAP-rifampin interactions and confer rifampin resistance in leprosy. AUTHOR SUMMARY: The emergence of primary and secondary drug resistance to rifampin in leprosy is a growing concern and poses a threat to the leprosy control and elimination measures globally. In the absence of an effective in-vitro system to detect and monitor phenotypic resistance to rifampin in leprosy, diagnosis mainly relies on the presence of mutations in drug resistance determining regions of the rpoB gene that encodes the ß-subunit of RNAP in M. leprae. Few labs in the world perform mouse food pad propagation of M. leprae in the presence of drugs (rifampin) to determine growth patterns and confirm resistance, however the duration of these methods lasts from 8 to 12 months making them impractical for diagnosis. Understanding molecular mechanisms of drug resistance is vital to associating mutations to clinically detected drug resistance in leprosy. Here we propose an in-silico saturation mutagenesis approach to comprehensively elucidate the structural implications of any mutations that exist or that can arise in the ß-subunit of RNAP in M. leprae. Most of the predicted mutations may not occur in M. leprae due to fitness costs but the information thus generated by this approach help decipher the impacts of mutations across the structure and conversely enable identification of stable regions in the protein that are least impacted by mutations (mutation coolspots) which can be a potential choice for small molecule binding and structure guided drug discovery.

Heteroexpression of Mycobacterium leprae hypothetical protein ML0190 provides protection against DNA-alkylating agent methyl methanesulfonate.

Repair of DNA alkylation damage is essential for maintaining genome integrity and Fe(II)/2-oxoglutarate(2OG)-dependent dioxygenase family of enzymes play crucial role in repairing some of the alkylation damages. Alkylation repair protein-B (AlkB) of Escherichia coli belongs to Fe(II)/2OG-dependent dioxygenase family and carries out DNA dealkylation repair. We report here identification of a hypothetical Mycobacterium leprae protein (accession no. ML0190) from the genomic database and show that this 615-bp open reading frame encodes a protein with sequence and structural similarity to Fe(II)/2OG-dependent dioxygenase AlkB. We identified mRNA transcript of this gene in the M. leprae infected clinical skin biopsy samples isolated from the leprosy patients. Heterologous expression of ML0190 in methyl methane sulfonate (MMS) sensitive and DNA repair deficient strain of Saccharomyces cerevisiae and Escherichia coli resulted in resistance to alkylating agent MM. The results of the present study imply that Mycobacterium leprae ML0190 is involved in protecting the bacterial genome from DNA alkylation damage.

Identification of novel open reading frames in the intergenic regions of Mycobacterium leprae genome and detection of transcript by qRT-PCR.

Mycobacterium leprae is an unculturable obligate pathogen and causative agent for debilitating human disease leprosy. Due to reductive genome evolution M leprae genome harbours large number of pseudogenes and small number of genes (∼1600 genes and ∼1300 pseudogenes). How M leprae remained a successful human parasite with small set of genes remains poorly understood and provided us the impetus to investigate the intergenic regions of M leprae genome for the presence of possible open reading frames (ORFs). In this work, we have manually scanned all the intergenic regions of M leprae genome and identified 106 potential ORFs. Among these, 12 are large ORFs: encoding hypothetical proteins (HP) of more than 100 amino acids. We have also found 67 ORFs encoding 50-100 amino acids proteins and another 27 ORFs for 30-50 amino acids peptides. We have validated the presence of transcripts for large HPs by quantitative reverse transcriptase PCR (qRT-PCR). Our results suggest that some of the M leprae large HPs are indeed expressed at low level in leprosy patients. The present results will shed light on the intergenic ORFs of M leprae and further our understanding of the pathogenesis of leprosy.

Publisher Correction: Structural Implications of Mutations Conferring Rifampin Resistance in Mycobacterium leprae.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Structural Implications of Mutations Conferring Rifampin Resistance in Mycobacterium leprae.

The rpoB gene encodes the ß subunit of RNA polymerase holoenzyme in Mycobacterium leprae (M. leprae). Missense mutations in the rpoB gene were identified as etiological factors for rifampin resistance in leprosy. In the present study, we identified mutations corresponding to rifampin resistance in relapsed leprosy cases from three hospitals in southern India which treat leprosy patients. DNA was extracted from skin biopsies of 35 relapse/multidrug therapy non-respondent leprosy cases, and PCR was performed to amplify the 276 bp rifampin resistance-determining region of the rpoB gene. PCR products were sequenced, and mutations were identified in four out of the 35 cases at codon positions D441Y, D441V, S437L and H476R. The structural and functional effects of these mutations were assessed in the context of three-dimensional comparative models of wild-type and mutant M. leprae RNA polymerase holoenzyme (RNAP), based on the recently solved crystal structures of RNAP of Mycobacterium tuberculosis, containing a synthetic nucleic acid scaffold and rifampin. The resistance mutations were observed to alter the hydrogen-bonding and hydrophobic interactions of rifampin and the 5' ribonucleotide of the growing RNA transcript. This study demonstrates that rifampin-resistant strains of M. leprae among leprosy patients in southern India are likely to arise from mutations that affect the drug-binding site and stability of RNAP.

Sequence homology and expression profile of genes associated with DNA repair pathways in Mycobacterium leprae.

BACKGROUND: Survival of Mycobacterium leprae, the causative bacteria for leprosy, in the human host is dependent to an extent on the ways in which its genome integrity is retained. DNA repair mechanisms protect bacterial DNA from damage induced by various stress factors. The current study is aimed at understanding the sequence and functional annotation of DNA repair genes in M. leprae. METHODS: T he genome of M. leprae was annotated using sequence alignment tools to identify DNA repair genes that have homologs in Mycobacterium tuberculosis and Escherichia coli. A set of 96 genes known to be involved in DNA repair mechanisms in E. coli and Mycobacteriaceae were chosen as a reference. Among these, 61 were identified in M. leprae based on sequence similarity and domain architecture. The 61 were classified into 36 characterized gene products (59%), 11 hypothetical proteins (18%), and 14 pseudogenes (23%). All these genes have homologs in M. tuberculosis and 49 (80.32%) in E. coli. A set of 12 genes which are absent in E. coli were present in M. leprae and in Mycobacteriaceae. These 61 genes were further investigated for their expression profiles in the whole transcriptome microarray data of M. leprae which was obtained from the signal intensities of 60bp probes, tiling the entire genome with 10bp overlaps. RESULTS: It was noted that transcripts corresponding to all the 61 genes were identified in the transcriptome data with varying expression levels ranging from 0.18 to 2.47 fold (normalized with 16SrRNA). The mRNA expression levels of a representative set of seven genes ( four annotated and three hypothetical protein coding genes) were analyzed using quantitative Polymerase Chain Reaction (qPCR) assays with RNA extracted from skin biopsies of 10 newly diagnosed, untreated leprosy cases. It was noted that RNA expression levels were higher for genes involved in homologous recombination whereas the genes with a low level of expression are involved in the direct repair pathway. CONCLUSION: This study provided preliminary information on the potential DNA repair pathways that are extant in M. leprae and the associated genes.

Analysis of a novel multiplex polymerase chain reaction assay as a sensitive tool for the diagnosis of indeterminate and tuberculoid forms of leprosy.

OBJECTIVE/BACKGROUND: Clinical diagnosis of indeterminate and tuberculoid leprosy is often difficult due to limited and confounding signs and symptoms. In the current study, we evaluated the utility of new multiplex polymerase chain reaction (PCR) using Mycobacterium leprae-specific DNA sequences in the pseudogene regions of ML1545, ML2180, and ML2179 for PCR-based diagnosis of indeterminate leprosy (IND) and leprosy cases across the immunological spectrum. The sensitivity was compared with that of RLEP PCR. METHODS: DNA was extracted from paraffin-embedded skin biopsy specimens of 220 leprosy cases, which were divided into IND (41), tuberculoid form (3), borderline tuberculoid (42), midborderline (3), borderline lepromatous (n=59), and lepromatous leprosy (72) cases. PCR positivity of both multiplex and RLEP PCR were compared in all the samples. A decision tree was constructed using the classification and regression trees algorithm to predict the probability of PCR positivity with the new multiplex PCR scheme in various clinical groups of leprosy. Sensitivity of each pseudogene target was determined using real-time PCR assays, and specificity was confirmed by PCR amplification of DNA extracted from three other mycobacterial species and skin biopsies of 44 non-leprosy cases. RESULTS: A multiplex PCR positivity of 75.61% was noted in IND cases when compared to that of 58.54% using RLEP PCR (P < 0.05). Enhanced multiplex PCR positivity was noted across various clinical groups in comparison to RLEP PCR. The decision tree classifier has predicted statistically significant probability for multiplex PCR positivity among RLEP-PCR negative group and clinical groups with a low bacillary load. CONCLUSION: This new multiplex PCR scheme can support the diagnosis of indeterminate and tuberculoid forms of leprosy with limited clinical manifestations and can be implemented in basic clinical/diagnostic setting that possess conventional PCR facilities.

Genomic diversity in Mycobacterium leprae isolates from leprosy cases in South India.

OBJECTIVE: The Objective of this study was to identify the strain diversity of Mycobacterium leprae in terms of SNP types and subtypes stratified as per genomic single nucleotide polymorphisms, in clinical isolates of leprosy patients from a tertiary care leprosy center in South India. Further, the associations of SNP types with clinical outcomes in leprosy were also investigated. METHODS: DNA was extracted from excisional skin biopsies of a total of 172 newly diagnosed untreated leprosy patients from a clinic in Tamil Nadu, in south India, that also serves patients from neighboring states. All the leprosy patients were those who voluntarily reported at the clinic during the study period of one year i.e., 2015. Clinical and histopathological details were collected at diagnosis and leprosy was confirmed through bacteriological smear examination and PCR for M. leprae specific RLEP region. SNP types and subtypes were determined by PCR amplification and Sanger sequencing of PCR products. RESULTS: M. leprae specific RLEP gene amplification was achieved in 160 out of 172 patients. Among 160 specimens 118(73.75%) were type 1 and 42 (26.25%) were type 2 and on subtyping it was noted that 88/160 (55.00%) were 1D, 25/160 (15.62%) 1C, 5/160 (3.12%) 1A, 33/160 (20.62%) 2G and 9/160 (5.62%) were 2H. CONCLUSION: Our results indicated that subtype 1D is predominant in the south Indian population. We also noted 2G, 1C and 1A in the patient sample tested. Additionally we identified subtype 2H for the first time in India.

Mycobacterium leprae specific genomic target in the promoter region of probable 4-alpha-glucanotransferase (ML1545) gene with potential sensitivity for polymerase chain reaction based diagnosis of leprosy.

With the absence of an effective diagnostic tool for leprosy, cases with negative bacteriological index and limited clinical manifestations often pose diagnostic challenges. In this study, we investigated the utility of a novel Mycobacterium leprae specific 112-bp DNA sequence in the promoter region of probable 4-alpha-glucanotransferase (pseudogene, ML1545) for polymerase chain reaction (PCR) based diagnosis of leprosy in comparison to that of the RLEP gene. DNA was extracted from slit skin scrapings of 180 newly diagnosed untreated leprosy cases that were classified as per Ridley Jopling classifications and bacteriological index (BI). Primers were designed using Primer Blast 3.0 and PCR was performed with annealing temperatures of 61°C for ML1545 and 58°C for the RLEP gene using conventional gradient PCR. The results indicated a significant increase in PCR positivity of ML1545 when compared to RLEP across the study groups (164/180 [91.11%] were positive for ML1545 whereas 114/180 (63.33%) were positive for RLEP [p<.0001, z=6.3]). Among 58 leprosy cases with negative BI, 28 (48.28%) were positive for RLEP and 48 (82.76%) were positive for ML1545 (p=.0001, z=3.8). Of the 42 borderline tuberculoid leprosy cases, 23 (54.76%) were positive for RLEP whereas 37 (88.09%) were positive for ML1545 (p<.0001, z=3.9). Increase in PCR positivity for ML1545 was also noted in lepromatous leprosy and BI-positive groups. ML1545 can be a potential gene target for PCR-based diagnosis of leprosy especially in cases where clinical manifestations were minimal.

Computational Modelling of Dapsone Interaction With Dihydropteroate Synthase in Mycobacterium leprae; Insights Into Molecular Basis of Dapsone Resistance in Leprosy.

The molecular basis for determination of resistance to anti-leprosy drugs is the presence of point mutations within the genes of Mycobacterium leprae (M. leprae) that encode active drug targets. The downstream structural and functional implications of these point mutations on drug targets were scarcely studied. In this study, we utilized computational tools to develop native and mutant protein models for 5 point mutations at codon positions 53 and 55 in 6-hydroxymethyl-7, 8-dihydropteroate synthase (DHPS) of M. leprae, an active target for dapsone encoded by folp1 gene, that confer resistance to dapsone. Molecular docking was performed to identify variations in dapsone interaction with mutant DHPS in terms of hydrogen bonding, hydrophobic interactions, and energy changes. Schrodinger Suite 2014-3 was used to build homology models and in performing molecular docking. An increase in volume of the binding cavities of mutant structures was noted when compared to native form indicating a weakening in interaction (60.7 Å(3) in native vs. 233.6 Å(3) in Thr53Ala, 659.9 Å(3) in Thr53Ile, 400 Å(3) for Thr53Val, 385 Å(3) for Pro55Arg, and 210 Å(3) for Pro55Leu). This was also reflected by changes in hydrogen bonds and decrease in hydrophobic interactions in the mutant models. The total binding energy (ΔG) decreased significantly in mutant forms when compared to the native form (-51.92 Kcal/mol for native vs. -35.64, -35.24, -46.47, -47.69, and -41.36 Kcal/mol for mutations Thr53Ala, Thr53Ile, Thr53Val, Pro55Arg, and Pro55Leu, respectively. In brief, this analysis provided structural and mechanistic insights to the degree of dapsone resistance contributed by each of these DHPS mutants in leprosy.